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ABSTRACT 


This  report  considers  the  peak  response  behavior  of  a  mass 
excited,  linear  mechanical  oscillator  when  the  applied  random 
excitation  is  either  stationary  or  nonstationary.  The  stationary 
random  excitation  is  Gaussian  bandwidth  limited  white  noise  and 
the  nonstationary  random  excitation  is  Gaussian  bandwidth  limited 
white  noise  shaped  in  time  by  (l)a  rectangular  envelope  function  or 
(2)  a  half-sine  envelope  function.  Thus,  the  nonstationary  excitation 
appears  as  a  pulse  of  white  noise  shaped  as  either  a  rectangle  or  a 
half- sine. 

Using  data  from  an  analog  computer  study,  two  topics  are 
explored  in  detail, 

•  the  expected  maximum  peak  response  in  a  finite  time 
interval  for  a  specified  probability  of  occurrence 

•  a  measure  of  the  time  duration  for  the  oscillator  to 
achieve  stationarity  in  its  response 

The  peak  response  is  expressed  as  parametric,  dimensionless  plots 
of  the  peak  to  rms  response  ratio  and  a  dimensionless  time  param¬ 
eter.  In  this  form,  the  oscillator  response  to  stationary  white  noise 
and  to  the  pulsed  random  excitation  can  be  conveniently  compared 
and  canbe  used  to  estimate  directly  the  time  duration  for  an  oscil¬ 
lator  to  achieve  stationarity  in  its  response.  It  is  found  that  for 
values  of  the  dimensionless  time  parameter  greater  than  about  one, 
the  stationary  results  can  be  applied  to  conservatively  predict  the 
oscillator  peak  response  to  the  pulsed  excitation. 

For  both  stationary  and  nonstationary  excitations,  the  peak 
behavior  is  noted  to  be  dependent  on  the  oscillator  damping  for  a 
small  number  of  response  cycles.  As  the  number  of  response 
cycles  becomes  large,  the  peak  to  rms  response  behavior  tends  to 
become  independent  of  damping.  v 
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THE  MECHANICAL  RESPONSE  OF  A 
LINEAR  MECHANICAL  OSCILLATOR  TO  STATIONARY 
AND  NONSTATIONARY  RANDOM  EXCITATION 


1.  INTRODUCTION 

This  report  treats  the  response  behavior  of  a  linear,  time  invariant 
mechanical  oscillator  subjected  to  Gaussian  random  excitation.  For  these 
conditions,  the  response  is  random  in  nature  and  is  characterized  by  Gaussian 
properties.  As  such,  only  the  first  two  statistical  moments  are  required  to 
provide  the  probabilistic  description  of  the  random  response.  The  complexity 
of  the  statistical  moments,  however,  depends  upon  whether  the  random  excita¬ 
tion  is  stationary  or  nonstationary,  whether  the  system  is  unimodal  or  multi¬ 
modal,  and  whether  the  transient  response  of  the  system  is  included  or  ignored. 

The  classic  unimodal  system  is  the  mechanical  oscillator  (also  commonly 
referred  to  as  a  single  degree-of-freedom  mechanical  system)  whereas  the 
multimodal  system  implies  a  distributed  elastic  structure.  The  distributed 
structure  generally  requires  attention  be  given  to,  space-time  correlation  of 
the  random  loading  whereas  the  mechanical  oscillator  requires  only  correlation 
in  time.  Including  transient  conditions  in  a  solution  introduces  nonstationarity 
into  the  response  for  the  initial  intervals  of  time. 

Reports  dealing  with  response  properties  categorically  include  the  mean 
value,  mean  square  response,  correlation  functions  and  power  spectra.  This 
report,  however,  comments  only  in  passing  on  these  properties  and  emphasizes 
the  peak  response  behavior  of  the  oscillator  when  the  applied  excitation  is 
(1)  stationary  and  (Z)  nonstationary. 

The  output  response  to  stationary  Gaussian  random  excitation  of  a  lightly 
damped,  linear  mechanical  oscillator  appears  as  sketched  in  Figure  1.  The 
output  is  Gaussian  and  appears  like  a  sinusoid  at  frequency  f  with  a  slowly 
varying  random  amplitude  and  random  phase.  The  amplitude  of  the  response 
is  shown  as  enveloped  in  time  so  that  it  becomes  proper  to  speak  of  the  statistical 
properties  of  the  envelope  as  well  as  the  response  time  history.  The 
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amplitude 


Figure  1.  Response  of  a  Lightly  Damped  Mechanical  Oscillator 
to  Broadband  Stationary  Random  Excitation 


instantaneous  probability  density  for  the  amplitude  of  the  response  time  history 
is  Gaussian.  The  probability  density  is  nearly  Rayleigh  for  the  distribution 
of  the  peaks  of  the  response  time  history.  As  the  oscillator  damping  approaches 
zero,  the  probability  density  of  the  peaks  of  the  time  history  approaches  a 
Rayleigh  distribution.  In  the  limit,  that  is,  when  the  oscillator  damping 
becomes  infinitesimally  small  thus  defining  an  infinitely  narrow  bandpass  filter, 
the  peaks  of  the  response  time  history  are  distributed  according  to  the  Rayleigh 
distribution.  The  probability  density  for  the  envelope  maxima  obtained  by 
passing  white  noise  through  an  ideal  bandpass  filter  is  given  by  Rice  on  page  223 
of  Reference  1.  Although  resembling  the  Rayleigh  distribution,  the  envelope 
maxima  density  is  noted  to  appear  as  neither  Rayleigh  nor  Gaussian.  Even 
though  these  properties  and  other  such  relationships  (Reference  2)  are  known 
for  simple  linear  systems  and  stationary  random  excitation,  several  important 
questions  still  cannot  be  answered. 

One  such  question  is,  "What  is  the  probability  of  exceeding  a  stated 
amplitude  within  a  finite  period  of  time?M  Expresses  in  an  alternate  manner, 
"For  a  stated  probability,  how  large  an  amplitude  will  occur  (on  the  average) 


2 


within  a  finite  period  of  time?"  To  answer  this  question  demands  a  time 
dependent  probability  whose  analytical  expression  still  remains  unknown  at 
this  writing.  The  solution,  however,  finds  immediate  application  in  the 
mechanical  design  of  missiles  and  spacecraft.  Since  exposure  to  maximum 
levels  of  random  excitation  occurs  for  short  time  periods  (for  example,  at 
launch,  near  Mach  1,  or  at  maximum  Q  in  the  flight  profile),  it  may  be 
possible  to  save  considerable  weight  by  establishing  structural  reliability 
criteria  based  on  a  time  related  probability  of  extreme  loads. 

Another  related  question  is,  "What  is  the  time  duration  for  a  mechanical 
oscillator  to  achieve  stationarity  in  its  response  to  impulsively  applied  random 
excitation?"  If  the  nonstationary  excitation  is  white  noise  enveloped  in  time 
by  the  unit  step  function,  the  response  will  not  be  stationary  for  the  initial 
intervals  of  time.  After  some  finite  time  duration,  however,  the  oscillator 
response  exhibits  stationary  characteristics  and  the  results  from  stationary 
analyses  then  can  be  applied.  If  the  nonstationary  excitation  is  white  noise 
enveloped  in  time  by  a  rectangular  or  half-sine  pulse,  the  response  charac¬ 
teristics  become  somewhat  more  difficult  to  codify.  A  solution  for  this  problem 
finds  application  in  the  design  of  systems  subjected  to  random  shock  loadings. 

Answers  to  these  questions  based  on  the  peak  response  of  a  linear 
oscillator  require  knowledge  of  the  distribution  of  the  response  maxima  in  time. 
Specifically,  it  is  required  to  know  the  distribution  of  time  between  level  cross¬ 
ings  or,  alternatively,  the  distribution  of  time  intervals  between  peaks  above  a 
stated  level.  This  is  an  extremely  difficult  problem  and  no  general  analytical 
solution  is  available  at  this  time  for  these  time  dependent  distributions.  Thus, 
it  becomes  appropriate  to  make  rather  sweeping  simplifying  assumptions  so  that 
an  approximate  time  dependent  distribution  can  be  written  for  the  response 
maxima.  The  main  portion  of  this  report  traces  the  analytical  development  of 
one  such  "approximate"  time  dependent  distribution  and  compares  predicted 
analytical  results  with  empirical  data  from  an  analog  computer  study. 
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2.  STATISTICAL  PROPERTIES 


Although  the  statistical  properties  of  a  Gaussian  random  process  are 
discussed  elsewhere,  it  is  instructive  to  briefly  consider  here  the  definitions 
of  those  functions  which  are  alluded  to  in  this  report.  The  statistical  quantities 
of  interest  are  the  mean,  the  variance,  and  the  covariance  functions.  These 
are  defined  for  stationary  and  nonstationary  random  processes. 

2.  1  STATIONARY  AND  ERGODIC  CONDITIONS 

Consider  the  statistical  moments  of  a  stationary  ergodic  random  process 
x(t).  The  assumption  of  ergodicity  allows  one  to  calculate  the  statistical 
properties  of  the  random  process  by  time  averaging  a  single  record  from  the 
random  process.  The  first  order  statistical  moment  is  the  mean  p  and  is 
defined  as 

T 

p  =  E  [  x(t)]  =  x  =  lim  ^  f  x(t)  dt  (2.1) 

T-*-oo  J0 


A  second  order  statistical  moment  is  the  autocorrelation  function 
and  is  defined  as 


R  (t) 

XX 


R  (t)  =  E  [  x(t)x(t  +  x)]  (2.2) 

2 

For  t  =  0,  the  autocorrelation  function  reduces  to  the  mean  square  value  ^ 
which  may  be  written  as 


R  (0)  =  iJj2  =  E[  x2(t)  ]  =  x2 

XX 


T 

lim  j  x2(t)  dt 

T-*oo  J0 


(2.3) 


The  covariance  function  C  (t)  also  is  a  second  order  statistical  moment  and 

xx 


appears  as 


For  r=  0,  the  covariance  function  reduces  to  the  variance  <r  which  may  be 
written  as 


C  (0)  =  o 
xx 


=  E 


(Ut)  -  [i] 


[x(t)  -  p]  2 


T 

lim  i  f  [x(t)  -  |jl]  2  dt  (2.5) 
T— ►oo  J0 


For  a  zero  mean,  the  covariance  function  is  identical  to  the  autocorrelation 
function,  i.  e.  , 


R  (t)  =  C  (t)  when  u  -  0  (2,6) 

xx  xx 


The  cross -covariance  function 
and  is  defined  as 


C  (r) 
xy 


is  still  another  second  order  moment 


C 

xy 


(t)  =  E 


c 


[x(t)  -  \l  ]  [ 

X 


y(t)  +  t)  - 


(2.7) 


For  zero  mean  values,  Eq. 
which  appears  as 


(2.7)  reduces  to  the  cross -correlation  function 


R 

xy 


(t)  =  E 


f 


x(t)]  [y(t  +  T)] 
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(2.8) 


Consequently,  for  zero  mean  values,  the  cross -covariance  function  and  the 

cross-correlation  function  are  identities.  R  (0)  has  no  particular  significant 

xy 

statistical  meaning. 

The  frequency  composition  of  x(t)  may  be  described  by  the  Fourier 
transform  of  the  autocorrelation  function  which  appears  as 


oo 


S  (f)  = 

XX  J  XX 

-oo 


f  .  .  .  -i2irfT 

I  R._. _(T)  e  dT 


(2.9) 


where  S  (f)  is  shown  in  terms  of  cyclical  frequencies  and  is  called  the  power 


xx 


spectral  density  function  with  the  units  of  mean  square  value  per  cps.  As 
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defined  here,  S^^(f)  is  a  two-sided  function  where  f  ranges  from  -oo  to  oo, 
and  is  related  to  the  spectral  density  in  radians  per  second  as  S  (f)  =  2tt  S  (oo) 

XX  XX 

Since  the  autocorrelation  function  is  symmetric  with  t,  then  the  imaginary 
part  of  Eq.  (2.9)  may  be  deleted  so  that 

oo 

S  (f)  =  I  R  (t)  cos  27rfT  dT  (2.  10) 

XX  y  XX 

“OO 

The  cross-spectrum  may  be  written  as 

oo 

S  (f)  =  /  R  (t)  e‘i2lTfT  dr  (2.11) 

xy  J  xy  '  7 

-oo  7 

The  real  part  of  the  cros s -spectrum  is  called  the  co-spectrum  and  the 
imaginary  part  of  the  cross  -  spectrum  is  called  the  quad-spectrum. 

Given  that  y(t)  defines  the  output  response  of  a  linear  time  invariant 
mechanical  system  due  to  the  input  excitation  x(t),  the  relationship  between 
the  input  and  response  power  spectral  densities  is  given  by 

Syy(f)  =  | H(f)|2  Sxx(f)  (2.  12) 

where  H(f)  is  the  frequency  response  function  for  the  linear  system.  Similarly, 
the  cross-spectrum  is  related  to  the  excitation  power  spectral  density  by 

Sxy(f)  =  H(f)  S^f)  (2.  13) 

Other  such  relationships  are  given  in  Reference  3. 

For  distributed  elastic  structures,  the  preceding  statistical  definitions 
must  be  modified  to  account  for  spatial  dependence  as  the  response  and 
excitation  are  both  functions  of  space  as  well  as  time.  Rather  than  dwell  at 
this  point  on  the  statistical  properties  of  distributed  structures,  the  reader 
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is  directed  to  Reference  4.  For  a  practical  interpretation  of  these  statistical 
properties  as  applied  to  the  testing  of  structures,  the  reader  is  directed  to 
References  5  and  6. 

2.2  NONSTATIONARY  CONDITIONS 

As  contrasted  to  the  statistical  moments  for  ergodic  random  processes, 
the  statistical  moments  for  nonstationary  processes  are  time  dependent  and 
are  calculated  by  averaging  over  an  ensemble  of  records  instead  of  time 
averaging  a  single  record.  Time  averaging,  however,  can  be  performed  but 
the  interpretation  of  these  results  is  not  obvious  as  is  discussed  in  References 
7  and  8. 

Assuming  that  x(t)  is  from  a  nonstationary  random  process,  the  mean 
is  given  as 


M-(t)  =  E[x(t)] 


The  autocorrelation  function  appears  as 

R  (t  ,t  )  =  E[x(t  )  x(t  )] 
xx  12  12 


(2.  14) 


(2.  15) 


Letting  t^  =  t^  =  t, 
value  defined  as 


the  autocorrelation  function  reduces  to  the  mean  square 


4>2(t)  =  E[x2(t)] 


(2.  16) 


The  covariance  function  may  be  written  as 


C 

xx 


{tl’t2) 


^xftj)  -  M-(t1)_]  [x(t. 


(2.  17) 


Letting  t^  =  =  t,  the  covariance  function  reduces  to  the  variance  of  x(t) 

which  may  be  written  as 


2 

cr 


(' 


(t)  =  E<[x(t)  -  p(t)] 


2] 


(2.  18) 
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For  a  zero  mean  value  over  all  time,  the  covariance  function  becomes  identical 
to  the  autocorrelation  function;  i.  e.,  when  |j,(t)  =  0, 


For  zero  mean  values,  the  cross-covariance  function  reduces  to  the  cro 

correlation  function  R  (t  ,  t  )  which  may  be  expressed  as 

xy  x  l 


V‘i  •  y  = E 


{[  x(ti 


)]  [y(t2>] 


(2.21) 


Basic  relationships  between  the  input  and  output  power  spectral  density,  the 
cross-spectral  density,  and  the  frequency  response  functions  for  a  linear 
system  are  discussed  in  Section  5  of  Reference  9. 


3.  RESPONSE  OF  A  MECHANICAL  OSCILLATOR 


This  section  deals  with  approximate  solutions  to  the  two  questions  cited 
in  the  introduction  regarding  the  oscillator  response  to  stationary  white  noise 
and  pulsed  random  excitation.  It  is  recalled  that  the  first  question  deals  with 
the  maximum  response  attained  in  time  T  with  probability  PM(P)  whereas 
the  second  question  considers  the  time  duration  for  the  oscillator  to  effectively 
attain  stationarity  in  its  response. 

As  one  answer  to  the  second  question,  Section  3.  1  paraphrases  the  work 
of  Caughey  and  Stumpf  (Reference  10)  where  the  excitation  is  white  noise  shaped 
by  the  unit  step  function.  Section  3.  2  considers  an  analog  study  by  Bar  noski 
and  MacNeal  (Reference  11)  which  provides  an  empirical  solution  to  the  first 
question  when  the  excitation  is  stationary.  Section  3.  3  discusses  analog  data 
by  Barnoski  and  MacNeal  (Reference  12)  which  provide  empirical  solutions  to 
both  questions  for  nonstationary  random  excitation. 

3.  1  RESPONSE  TO  WHITE  NOISE  SHAPED  IN  TIME  BY  THE  STEP  FUNCTION 

Much  of  the  literature  concerning  the  response  of  simple  mechanical 
systems  to  stationary  random  excitation  considers  only  the  stationary  aspects 
of  the  dynamic  response.  In  including  the  transient  motion  of  the  mechanical 
oscillator,  the  response  of  the  system  to  white  noise  shaped  by  the  step  function 
is  nonstationary  for  the  initial  intervals  of  time.  Two  excellent  papers  in  this 
subject  area  are  found  in  References  10  and  13. 

In  keeping  with  the  discussion  of  Reference  10,  the  equation  of  motion 
for  the  mechanical  oscillator  acted  upon  by  a  random  acceleration  excitation 
applied  to  the  mass  may  be  written  as 

x  +  2c o  £x  +  co  ^x  =  a(t)  (3.  1) 

n  n 


where 
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t,  -  damping  ratio  of  the  mechanical  oscillator 
x  =  displacement  of  the  mass  from  static  equilibrium 

x  =  velocity  of  the  mass 

•  # 

x  =  acceleration  of  the  mass 

wn  ~  undamped  natural  frequency  of  the  mechanical  oscillator 
=  random  acceleration  acting  on  the  mass 

The  undamped  natural  frequency  and  damping  ratio  are  related  to  the  physical 
parameters  of  the  mechanical  system  by 


2 

co 


n 


_k 

m 


(3.2) 


6  = 


c 


2m  co 

n 


where 

k  =  linear  spring  constant 

m  =  mass  of  the  system 

c  =  viscous  damping  coefficient 

The  random  excitation  is  assumed  to  have  the  following  properties 

•  a(t)  is  stationary 

•  a(t)  is  Gaussian 

•  ot(t)  has  a  zero  mean  value 

•  Qf(t)  has  the  power  spectrum  G^(u>)  which  is  considered  to 
be  smooth  and  contains  no  sharp  peaks 

•  <*(t)  is  assumed  to  be  mean  square  continuous 


(3.3) 
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Assuming  small  values  of  damping  (£,  <  1)  and  that  the  initial  conditions  are 


x(0)  =  a;  x(0)  =  b 


(3.4) 


then,  a  solution  of  (3.  1)  may  be  written  as 


x(t)  =  a  e 


-£co  t 
n 


&co 


cos  co-t  +  — 
a  a 


n  •  tl 

-  sm  co  t 


h  * 

.  o  n  . 

+  -  e  sin  co  ,t  + 

“d  d  -0 


f 

J  h(t  -  T )  a 


(t)  dT 


(3.5) 


where 


“d  -V 


1  -  co  =  damped  natural  frequency  of  the  oscillator  (3.  6) 


and  for  t  >  0 


-£co  t 
n 


h(t)  = 


co 


sin  go  t  —  impulse  response  of  the  system 
d 


(3.7) 


Since  the  response  is  noted  to  be  a  Gaussian  process,  only  the  mean  and  the 
variance  need  be  calculated  to  characterize  the  response  properties.  The  mean 
value  of  x(t)  is  formed  by  an  ensemble  average  and  appears  as 


-co  £t  " 

p(t)  =  E[x(t)]  =  a  e  n  < 


co  t 

cos  co,t  +  — —  sin  co^t 
d  d  d 


h  -co  Ct  d 

+  -  e  n  sin  w  f  +  I  h(t  -  t)  E[q?(t)]  dT  (3.  8) 

go  ,  d 


For  the  mean  value  of  the  excitation  asumed  equal  to  zero,  the  mean  of  the 


response  appears  as 


(3.  9) 


M-(t)  =  E[x(t)]  =  a  e 


-w  it 


< 


conC 

cos  ca  +  -  sin  go  t 

d  go,  d 


-u  £t 

,  b  n  . 

+  -  e  sin  go  t 

go  ,  d 

d 


The  variance  of  x(t)  is  formed  also  by  an  ensemble  average  and  is  given 
as 


2 

a* 


(t)  =  E 


«[ 


[x(t)  -  Ji(t)] 


5 


(3.  10) 


In  somewhat  different  form,  the  variance  for  the  mechanical  oscillator 
appears  as 


crZ(t)  = 


(L 


h(t  -  t)  h(t  -  t1  )  E  [  ^(t)  a(T’ )]  dT  dTf 


(3.  11) 


From  Eq.  (2.  2),  the  expression  containing  the  product  of  the  acceleration 
terms  is  noted  to  be  the  autocorrelation  function  for  cv(r).  This  function  may 
be  written  as 

R  (t ,  t  1 )  =  e[<*(t)  •  or(T')]  (3.  12) 

aa 

Note  that  the  t  nomenclature  in  Eq.  (2.2)  refers  to  a  time  difference  whereas 
the  t  symbols  in  Eq.  (3.11)  refer  to  variables  of  integration.  For  the 
stationary  assumption,  the  autocorrelation  function  depends  only  on  the  time 
difference  (t-t1)  rather  than  on  t  and  t1  individually.  Expressed  in  terms 
of  the  physically  realizable  one-sided  power  spectrum,  Eq.  (3.12)  appears 
as 

oo 

RaJAt)  =  RaJT  ~  T')  =  f  Ga(w)  cos  “(T  "  T')  dw  (3.13) 

Substituting  (3.13)  and  (3.7)  into  (3.11)  and  performing  the  double  time 
integration  yields 
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<r2(t) 


(co )  r  -2co^£,tr 

=  i  h^f2V  +  e  “  L 


1  +  -  L  sm  co  t  cos  co 

co  _  d  d 

d 


/  2^  \ 

e  (2  cos  co  -t  4*  -  sin  co^t]  cos  cot 

'  d  “d  d 


(3.  14) 


where 


2  2  2 

co  £t  _  (go  £)  -  co  -  +  co  ~ 

n  2co  .  .  ^  ,  n  d  .2j 

—  sm  co  nt  sm  cot  +  - - -  sm  co  -t 

co  .  d  2  d 

d  co  _ 

d 


|  Z  ( co)  |  =  ^co^  -  co  ^  +  (2  co  co^f,) 


(3.  15) 


If  G  ( co)  is  a  smooth  function  of  co  with  no  sharp  peaks  and  £  <<  1,  then 
a 

Eq.  (3.  14)  can  be  approximated  by  Laplace's  method  of  evaluating  integrals 
which  provides 


ttG  (  co  ) 

r2,t)~  -JL» 

4t,  co 

n 


-2co  It  r  n  n2 

e  n  2  (2con^  2 

1  -  - - -  co  n  +  - - -  sin  co  ,t+  co  co  ,tsin  ,t 

2  d  2  d  n  d  d 

CO 

d  u 


(3.  16) 


where  G  (co  )  denotes  the  magnitude  of  G  (co)  at  co  .  Alternatively, 
an  an 

Eq.  (3.  16)  can  be  expressed  as  a  nondimensional  variance  in  the  form 


2crZ(t)  “n  „  1 

TT  G  (co  )  2Ci 
a  n 


f  "^"n*  r  2|2  2  /  2\  2 

\  <■  i -i  \  /  a 


(-c2) 


sin  2  1  - 


f  vi 


(3.  17) 


Figure  2  is  a  family  of  curves  of  Eq.  (3.  17)  for  various  values  of  damping  £ 
and  is  an  extension  of  the  same  plot  shown  in  Reference  10.  All  curves 
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igure  2.  Dimensionless  Variance  of  the  Displacement  Response  of  a  Mechanical  Oscillator 
Excited  by  White  Noise  Shaped  in  Time  by  the  Step  Function 


begin  at  zero  and,  except  for  £,  =  0,  asymptotically  approach  1/2 £,  as  limiting 
values.  Of  particular  interest  is  the  value  of  co  t  as  the  magnitude  of  the 
ordinate  approaches  an  asymptotic  value.  This  value  yields  the  time  or 
number  of  cycles  for  the  mechanical  oscillator  to  attain  stationarity  in  its 
response.  For  example,  co^tA<#  6tt  when  £  =  0.10;  this  corresponds  to 
approximately  three  cycles  for  the  system  response  to  achieve  stationarity. 

3.  2  PEAK  RESPONSE  PROPERTIES  TO 
STATIONARY  RANDOM  EXCITATION 

The  probability  density  is  not  yet  known  for  the  distribution  in  time  of 
the  peak  response  of  a  mechanical  oscillator  excited  by  white  noise.  Rice 
(page  223  of  Reference  1),  however,  derived  a  probability  density  for  the 
envelope  maxima  of  the  output  from  an  ideal  rectangular  bandpass  filter 
excited  by  white  noise.  In  addition.  Rice  showed  that  the  average  number  of 
envelope  maxima  is 

H  =  0.641B  (3.18) 

n 

where  B  is  the  noise  bandwidth  of  the  filter,  Aspinwall  (Reference  14),  in 
turn,  defined  a  mechanical  equivalent  to  the  electrical  ideal  bandpass  filter 
and  used  these  two  results  of  Rice  to  obtain  an  approximate  answer  to  the 
question  regarding  the  probability  of  exceeding  a  stated  amplitude  within  a 
finite  time  interval.  For  completeness  of  this  discussion,  it  is  appropriate 
to  paraphrase  the  arguments  of  Aspinwall, 

An  equivalent  ideal  bandpass  filter  is  obtained  in  the  following  manner. 

In  Figure  3,  the  frequency  response  characteristics  are  shown  for  both  a 
rectangular  bandpass  filter  and  a  mechanical  oscillator.  The  ordinate  is  noted 
as  an  absolute  value  of  the  frequency  response  function  and  the  x-axis  is  the 
ratio  of  the  excitation  frequency  to  the  undamped  natural  frequency  of  the 
oscillator.  B  is  shown  as  the  bandwidth  for  an  ideal  rectangular  filter. 
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B 


Figure  3.  Magnitude  of  the  Frequency  Response  Function  for  a 
Mechanical  Oscillator  and  an  Ideal  Rectangular  Band¬ 
pass  Filter 

The  sketch  for  the  mechanical  oscillator  depicts  the  magnitude  of  the 
velocity  to  force  frequency  response  function  for  a  mass  excited  oscillator. 
This  function  is  expressed  as 


H(co)  - 


f(t) 


m 


oo 


s  +  Zt,  oo  + 

n  s 


n 


1/c 


1  + 


» 

m 

—  s 
c 


1  + 


ms 


(3.  19) 


where  s  -  iw  and  other  relationships  are  provided  by  Eqs.  (3.  2)  and  (3.  3). 
In  terms  of  the  frequency  in  cycles  per  second,  the  absolute  magnitude  of 
(3.  19)  is 


H(f)  = 


2ir  f/k 


V  hr]2  4  -01 


(3.  20) 
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As  shown  by  the  sketch  in  Figure  3,  Eq,  (3,20)  is  zero  at  both  extremes  of  the 
frequency  axis  and,for  small  values  of  damping,  has  a  maximum  value  given 
approximately  by 


I  H(f)l 


max 


TTf 

_ n 

u 


(3.21) 


Bandwidth  equivalence  between  the  ideal  bandpass  filter  and  the  mechani¬ 
cal  oscillator  is  defined  by  requiring  the  output  responses  to  white  noise  to 
have  the  same  mean  square  value.  It  is  tacitly  assumed  and  shown  in 
Figure  3  that  the  center  frequencies  and  gain  for  both  systems  are  the  same. 
Thus,  the  mean  square  responses  to  white  noise  are 


mechanical 

oscillator 


I  H(f )  1 


2 

max 


G„<f) 


(3.22) 


ideal 

bandpass 

filter 


x2  =  B  ♦  I  H(f)  I  2  •  G.(f) 

max  '  0 


(3.23) 


where  B  is  the  noise  bandwidth,  H(f)|  is  the  square  of  the  maximum 

max 

value  of  the  frequency  response  function,  and  G^( f)  is  the  magnitude  of  the 
white  noise  input  excitation.  Equating  (3.22)  and  (3.23)  yields  the  following 
relationship  between  the  filter  noise  bandwidth  and  the  oscillator  half-power 
point  bandwidth 


B  = 


(3.24) 


where  the  half-power  point  bandwidth  is 


(3.25) 


Substituting  Eq.  (3.24)  into  Eq.  (3.18)  yields  the  average  number  of  envelope 
maxima  as 


n 


H  ~  -rr 
n  —  Q 


(3.26) 
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Alternatively,  substituting  (3.21)  into  (3.22)  yields  the  mean  square  velocity 
response  as 

—  co 

•  2  Q  n 

X  =  4G0(£)~  <3-27> 

k 

Having  an  estimate  for  an  expression  for  predicting  the  probability 

of  exceeding  a  given  amplitude  within  a  finite  time  interval  is  argued  as 
follows.  If  it  is  assumed  that  the  envelope  maxima  are  independent  random 
variables  and  that  the  number  of  envelope  maxima  for  the  time  interval  T 
is  given  by  p^T  ,  THEN  the  occurrence  of  each  envelope  maxima  is  one  of 
M-nT  Bernoulli  trials  where  success  is  attained  if  the  maxima  <  (3  .  The 
probability  that  the  maximum  value  of  the  envelope  maxima  remains  <(3 
during  the  time  interval  T  is  then  given  as 

pM(P)  =  [pE(P>rn  (3.28) 

where 

Pr°bability  that  M,  the  maximum  value  for  the  envelope 
peaks  in  time  T,  is  (3  times  the  rms  value  of  the  response 

PE(P)  the  probability  that  the  maximum  value  of  any  envelope  peak 

is  times  the  rms  value  of  the  output  response.  For  the 
ideal  bandpass  filter,  this  value  is  obtained  by  integrating  the 
probability  density  function  of  the  envelope  maxima  (Reference  1, 
page  223)  over  the  limits  of  interest. 

^n^  *s  average  number  of  envelope  maxima  during  time  T 

Pn  is  the  average  number  of  envelope  maxima  per  second 

T  is  the  time  interval  for  observing  the  output  response 

Alternatively,  the  probability  of  an  envelope  maxima  exceeding  (3  in  time 
T  is  expressed  as 

V>P>  =  1  -  PM(P)  (3.29) 
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By  way  of  interpreting  Eqs.  (3.28)  and  (3.29)  consider  their  application 
to  this  classic  statistical  problem  of  flipping  a  single  coin;  what  is  the  proba¬ 
bility  of  flipping  five  heads  in  a  row?  For  this  problem,  Eq.  (3.28)  can  be 
used  directly.  In  this  context, 


pm»> 

PE(P) 


denotes  the  probability  of  flipping  five  heads  in  a  row 

denotes  the  probability  of  flipping  a  head  in  a  single 
toss  of  the  coin,  that  is  1/2 


denotes  the  total  number  of  trials  or  flips  of  the  coin, 
which,  in  this  example,  equals  five 


so  that 


PM(P)  =  (1/2)  =  0.031 


(3.30) 


From  Eq.  (3.29),  the  probability  of  not  flipping  five  heads  in  a  row  is 


simply 


P  (>p)  =  1  -  .  031  =  0.  969 


(3.31) 


Substituting  (3.26)  into  (3.28)  yields  P.^.((3)  f°r  the  mechanical  oscillator 


pmW  -  (PE«»] 


(3.32) 


where 


(3.33) 


T  is  noted  as  a  dimensionless  time  parameter  and  f^T  may  be  interpreted 
as  the  number  of  cycles  for  which  the  random  response  is  observed*  The 
P  (|3)  term  is  obtained  directly  from  the  probability  density  for  envelope 

hi 

maxima  given  by  Rice  on  page  223  of  Reference  1. 
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As  sketched  in  Figure  4,  Eq.  (3.32)  can  be  displayed  as  a  family  of 
probability  curves  in  P  (p)  plotted  as  log  t"  versus  p.  These  curves 
represent  the  results  of  a  theoretical  development  by  Aspinwall  and  are 
^PPropr iately  called  Aspinwall  curves  in  this  report.  It  is  remembered  that 
these  curves  are  applicable  to  the  peak  properties  of  the  stationary  response 
for  a  lightly  damped  mechanical  oscillator  where  the  excitation  is  stationary 
white  noise. 


P  - ► 

Figure  4.  Parametric  Plot  for  the  Peak  to  RMS  Response  of  a 
Mechanical  Oscillator  Subjected  to  Broadband 
Stationary  Random  Excitation 


As  mentioned  in  Reference  14,  the  errors  inherent  in  the  Aspinwall 
curves  are  associated  with  the  following  assumptions: 

(1)  Differences  between  the  actual  number  of  envelope  maxima 

and  the  approximate  number  given  as  \i  T. 

n 

(2)  Differences  in  p  for  the  rectangular  bandpass  filter  and 
the  mechanical  oscillator. 
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(3) 


Errors  in  assuming  the  P^((3)  developed  by  Rice  for  an  ideal 
rectangular  filter  to  be  directly  applicable  for  the  band¬ 
pass  of  the  mechanical  oscillator. 

(4)  Errors  in  assuming  the  envelope  maxima  to  be  statistically 
independent. 


Although  these  factors  bias  the  numerical  results  of  the  theoretical  analysis, 
the  trends  displayed  by  the  curves  in  Figure  4  were  believed  to  be  valid. 

In  Reference  11,  Barnoski  and  MacNeal  examined  the  validity  of 
Aspinwall  curves  bycomparing  these  theoretical  results  with  empirical  data 
from  an  analog  computer  study.  For  this  study,  a  mechanical  oscillator  was 
simulated  by  the  operational  amplifier  circuits  shown  as  Figure  5.  The 
triangles  denote  simple  electronic  amplifiers,  the  K's  represent  potentiometers, 
C!s  are  capacitors,  R’s  are  resistors,  E^  is  the  input  excitation  and  e 
represents  the  velocity  of  the  mechanical  system.  In  terms  of  the  s  operator, 
the  ratio  of  the  velocity  to  the  input  excitation  is 

KR 
3  3 


KZR4 


R3C1S 


R1C1R2C2S 


Comparing  the  terms  of  (3.19)  with  (3.34)  yields 


K3R3 

K2R4 


(3.35) 


R3C1 


(3.36) 


R1C1R2C2 


(3.37) 
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Selecting  the  appropriate  values  for  the  components  in  (3.35),  (3.36),  and 
(3.37)  guarantees  the  simulation  of  the  velocity  to  force  frequency  response 
function  for  the  mechanical  oscillator. 


Figure  5.  Operational  Amplifier  Circuits  for  a  Linear 
Mechanical  Oscillator 


In  the  analog  study,  the  oscillator  was  excited  continuously  by  white 

noise  and  the  velocity  response  was  observed  for  time  intervals  (T)  ranging 

from  a  minimum  of  0.1  second  to  a  maximum  of  20  seconds.  The  white  noise 

was  provided  by  a  random  noise  generator  and  the  noise  level  was  set  so 

^6cr  peaks  of  the  output  voltage  would  not  be  clipped  by  the  amplifiers. 

For  each  time  interval  and  value  of  Q  ,  the  output  signal  was  sampled 

100  times  where  the  time  duration  for  each  sampling  was  T  seconds. 

Figure  6  shows  two  typical  time  histories  as  observed  during  the  sampling 

interval  T=0.1  second,  and  where  Q  =  20  and  f  =  318  cps. 

n 

F or  each  sampling,  the  maximum  positive  and  negative  values  for  the 
output  signal  was  measured  and  recorded.  After  100  samplings,  these  data 
then  were  arranged  in  rank  order  to  yield  histograms  of  probability  versus  (3. 
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Through  each  histogram,  a  mean  curve  was  drawn  which  represents  an 
empirical  plot  of  P  (P)  versus  (3.  These  empirical  curves  were  used  to 
create  empirical  plots  of  f  T/Q  versus  (3  for  values  of  P^tP)  as  illustrated 
by  Figures  7  and  8. 

The  coordinates  in  these  figures  are  noted  to  be  the  coordinates  of 
Figure  4  so  that  the  Aspinwall  probability  curves  can  be  compared  directly 
with  the  empirical  curves.  The  empirical  curves  are  plotted  as  a  family  in 
Q  for  Q  values  of  5,  10,  20,  and  50.  Each  data  point  shown  was  obtained  from 
mean  curves  of  the  histograms  alluded  to  in  the  previous  paragraph.  Also 
shown  are  the  approximate  number  of  cycles  for  which  the  oscillator  was 
observed.  It  is  interesting  to  note  that  for  the  larger  values  of  f^T/Q,  the 
constant  probability  curves  tend  to  appear  as  straight  lines  on  the  semi¬ 
log  graph  paper.  This  implies  the  f^T/Q  and  (3  are  related  in  form  as 

f  T 

log  =  mp  _j_  constant  (3.38) 

where  m  is  the  slope  of  the  straight  lines. 

From  Figures  7  and  8,  it  can  be  concluded  that  the  empirical  data  lend 
support  to  the  theory  embodied  in  Aspinwall  curves  although  the  response 
ratios  ((3)  are  consistently  higher  than  those  predicted  by  the  Aspinwall  theory. 
In  other  words,  the  numerical  values  provided  by  the  theoretical  curves  are 
not  conservative.  By  noting  the  vertical  alignment  of  data  points  for  the 
larger  number  of  cycles,  it  may  be  concluded  that  the  response  ratio  ((3) 
tends  to  become  independent  of  Q  for  a  large  number  of  cycles. 

Figure  9  is  a  family  of  curves  in  P  (p)  for  tiie  Peak  to  rms  response 
ratio  ((3)  versus  the  dimensionless  time  parameter  (f^T/Q).  These  curves 
are  based  on  the  empirical  plots  from  Figures  7  and  8,  as  well  as  similar 
data  from  Reference  12.  Only  two  values  of  P  (P)  are  shown  and  each 
p  (P)  curve  is  noted  to  be  influenced  by  the  damping  values  of  the  oscillator. 
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The  p  ordinate  values  where  T  =  0  become  the  instantaneous  values  of  the 
amplitude  time  history.  Since  the  probability  density  of  the  instantaneous 
values  is  normal,  these  magnitudes  are  obtained  from  most  any  tabulated 
table  of  the  Gaussian  density  function.  All  the  curves  of  Figure  9  exhibit  a 
rapid  rise  in  the  p  value  for  initial  intervals  of  f  T  /  Q.  After  passing  through 
a  knee  in  the  vicinity  of  f^T/Q  =  1,  the  p  values  increase  very  slowly  for 
relatively  large  increases  in  f^T/Q.  Although  the  plots  are  shown  only  up 
to  f^T/ Q  =  7,  it  is  understood  that  these  curves,  in  the  limit,  extend  to 
infinity  where  the  p  ratio  theoretically  becomes  infinite.  Thus,  an  infinite 
peak  is  to  be  expected  in  an  infinite  time. 

The  curves  of  Figure  9  show  that  peak  to  rms  response  ratios  of 
3[pM(P)  =  °.6°]  and  4[PM(p)  =  0.95]  are  not  uncommon  in  short  time  inter¬ 
vals  for  a  linear  oscillator.  On  the  other  hand,  a  large  time  lapse  is 
required  before  p  ratios  of  5  or  6  are  expected.  These  curves  assume  signi¬ 
ficance  to  a  designer  if  they  are  applied  to  establish  structural  reliability 
criteria  based  on  a  time  dependent  probability  of  extreme  loads.  In  addition, 
these  curves  can  be  used  to  define  stationarity  in  the  peak  response  as  the 
shown  plots  are  the  stationary  results  of  an  oscillator  excited  by  stationary 
white  noise. 

Figure  10  depicts  families  of  curves  in  P  (P)  for  the  peak  to  rms 

response  ratio  plotted  versus  the  number  of  cycles  for  the  oscillator  (f  T). 

n 

The  equation  for  this  theoretical  boundary  was  developed  by  Thrall  using  non- 


parametric  statistics  (Reference  15)  and  can  be  expressed  as 


(T  [I-PM<P>]  P2/Z 

xi=  - - -  e 

n  Z 


(3.39) 


The  shown  theoretical  curve  is  only  for  P  (p)  =  0.95  whereas  the  remaining 
PM^  Plots  are  empirical  plots  formed  from  the  analog  data  of  Figures  7  and 
8.  Figure  10  clearly  points  out  that  the  (3  response  ratio  becomes  independent 


of  the  oscillator  damping  for  a  large  number  of  cycles,,  Although  the 
theoretical  P  ((3)  plot  is  based  upon  an  average  value  for  the  number  of 
level  crossings  per  unit  time,  Eq.  (3.39)  appears  as  a  valid  bound  for  large 
values  of  T.  For  small  T  values,  however,  the  theoretical  bound  is  noted  to 
be  inexact,  but  conservative,  when  compared  with  the  empirical  curves. 

By  comparing  the  Q  curves  for  any  value  of  ^  noticed  in 

Figure  10  that  the  (3  values  are  higher  for  the  larger  damping  values.  This 
behavior  requires  further  explanation.  The  output  maxima  response  of  an 
oscillator  excited  by  white  noise  may  be  considered  to  consist  of  clumps  of 
envelope  maxima  per  unit  time.  For  a  very  high  Q  system,  the  envelope 
maxima  (clumps)  are  noted  to  be  well  separated  in  time.  As  damping  is 
increased,  the  time  separation  between  envelope  maxima  becomes  less  so 
that  a  greater  number  of  maxima  occur  per  unit  time.  Hence,  the  greater 
the  clump  density  in  time,  the  greater  the  likelihood  of  experiencing  higher 
peak  to  rms  values. 
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Figure  6.  Two  Typical  Time  Histories  for  the  Response  of  a  Mechanical 
Oscillator  Excited  by  Stationary  White  Noise 
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Figure  7.  Dimensionless  Time  Parameter  Versus  the  Ratio  of  Peak  to 

RMS  Response  for  the  Output  of  a  Mechanical  Oscillator  Excited 
by  Stationary  White  Noise:  P^(P)  =  0.  60 
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Figure  8.  Dimensionless  Time  Parameter  Versus  Ratio  of  Peak  to  RMS 
Response  for  the  Output  of  a  Mechanical  Oscillator  Excited  by 

Stationary  White  Noise:  P  (P)  -  0-  95 
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Figure  9.  Experimental  Envelope  for  the  Peak  to  RMS  Response  of  the  Mechanical  Oscillator  Excited 
by  Stationary  White  Noise  ((3  versus  f  T/Q) 
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Figure  10.  Experimental  Results  for  the  Peak  to  RMS  Response  of  the  Mechanical  Oscillator  Excited 
by  Stationary  White  Noise  ((3  versus  f  t) 


3  ,  3  PEAK  RESPONSE  PROPERTIES  TO 

NONSTATIONARY  RANDOM  EXCITATION 

Attention  is  focused  in  this  section  upon  the  peak  response  behavior 
of  a  mechanical  oscillator  when  the  excitation  is  white  noise  shaped  in  time 
by  either  a  rectangular  or  half- sine  envelope  function.  Hence,  the  excitation 
is  nonstationary  and  appears  as  a  pulse  of  white  noise  shaped  either  as  a 
rectangle  or  a  half- sine.  This  problem  is  typically  associated  with  structural 
response  to  shock  loading  where  the  shock  can  be  due  to  phenomenon  as  earth¬ 
quakes,  nuclear  blasts,  or  missile  launch,  Barnoski  and  MacNeal  (Reference  12) 
treated  this  problem  for  a  linear  oscillator  using  an  analog  computer  and  much 
of  their  data  form  the  basis  of  this  discussion. 

To  display  the  peak  response  results  to  nonstationary  excitation  using 
the  f^T/Q  and  (3  parameters  of  the  previous  section,  a  time  duration 
analogous  to  the  sampling  interval  T  need  be  defined.  It  is  recalled  that 
T  refers  to  the  time  interval  for  sampling  the  response  of  the  linear  oscillator 
to  stationary  white  noise.  For  the  nonstationary  study,  the  response  is  con¬ 
tinuously  monitored  so  that  a  time  interval  analogous  to  T  is  stated  in  terms 
of  the  duration  of  the  input  pulse.  An  effective  or  analogous  time  interval  At 
is  defined  by  requiring  the  input  energy  in  the  nonstationary  pulse  to  be  equal 
to  an  equivalent  amount  of  energy  for  stationary  conditions.  Hence,  the  non¬ 
stationary  time  interval  is  expressed  as 

oo 

At  = - — -  I  E2(t)  dt  (3.40) 

E  J 0 

max 

where  E(t)  is  an  envelope  function  denoted  as 

E(t)  =  E  g(t)  (3,41) 

max 


31 


E  is  the  maximum  value  of  the  envelope  function  and  g(t)  is  the  shape 

max 

function.  For  the  rectangular  and  half- sine  shape  function,  E(t)  and  the 
effective  time  intervals  (At)  for  these  pulses  appear  as  shown  in  Figure  11. 


Rectangular 
Shape  Function 


Half-Sine 
Shape  Function 


E(t) 


time 


Figure  11.  Rectangular  and  Half-Sine  Envelope  Functions  for  the 
Pulsed  Random  Excitation, 


Thus,  the  effective  time  interval  for  the  rectangular  pulse  is  equal  to  its 
actual  time  duration  whereas  the  effective  time  interval  for  the  half-sine 
pulse  is  one -half  its  actual  time  duration.  It  is  to  be  noted  that  the  rectangular 
shape  function  approaches  a  step  function  as  At  becomes  very  large. 

The  dimensionless  time  parameter  for  the  nonstationary  conditions 
becomes 


where  Eq.  (3.42)  differs  from  Eq.  (3,33)  only  in  the  definition  of  the  time 
intervals  T  and  At. 
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The  expression  for  the  pulsed  excitation  can  be  written  as 


f(t)  =  G  E(t)  (3o43) 

where  is  the  magnitude  of  the  white  noise  spectrum.  Equation  (3,43)  is 

created  electrically  by  multiplying  the  outputs  from  a  random  noise  source 
and  an  envelope  generator. 

In  their  analog  study,  Barnoski  and  MacNeal  varied  T  in  discrete 
steps  ranging  from  a  minimum  value  of  0,  25  to  a  maximum  value  of  5.  00, 

This  maximum  value  of  T  is  considered  large  enough  so  that  At  is  long 
compared  to  the  natural  period  of  the  oscillator.  For  this  condition,  the 
input  excitation  appears  as  stationary  white  noise  to  the  oscillator  and  the 
(3  response  should  be  nearly  the  same  as  in  Figure  9.  The  undamped  natural 
frequency  of  the  system  (f  )  was  set  at  159  cps  and  Q  values  were  5,  20,  and 
50, 

sU 

For  a  specified  T  ,  the  oscillator  was  excited  by  either  a  rectangular 
pulse  or  a  half- sine  pulse  of  white  noise  and  the  maximum  positive  and 
negative  values  of  the  peak  to  rms  time  histories  were  recorded.  This  pro- 
cedure  was  repeated  100  times  for  each  T  value.  As  with  the  analog  data 
of  Section  3,  2,  these  recorded  data  were  arranged  in  rank  order  to  construct 
PM(P)  versus  (3  curves;  then,  to  create  the  desired  i^Ar/Q  versus  (3  plots. 
Figure  12  shows  a  typical  rectangular  pulse  of  white  noise  and  the 
resulting  response  of  a  linear  oscillator  to  such  an  excitation.  The  experi¬ 
mental  constants  are  Q=20  and  At  =  0.0314  second  which  defines 

T  =  f  At/Q=0,25,  The  effective  time  interval  (At)  defined  by  Eq.  (3,40)  is 
n 

equal  to  the  actual  time  duration  of  the  pulse  and  is  noted  to  be  five  times 
the  natural  period  of  the  oscillator.  The  response  is  a  velocity  time  history 
and  is  observed  as  a  damped  harmonic  oscillation  after  the  pulse  is  terminated. 

Figure  13  shows  a  half-sine  envelope  function  and  a  typical  resulting 
pulsed  random  excitation  when  this  envelope  function  is  multiplied  with  the 
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output  of  a  random  noise  generator.  This  half- sine  pulse  of  white  noise  has 

an  effective  time  interval  (At)  of  one -half  the  time  duration  of  the  envelope 

function.  As  with  the  rectangular  pulse,  the  effective  time  duration  is  five 

times  the  period  of  the  mechanical  oscillator.  Figure  14  depicts  two  typical 

responses  of  the  oscillator  when  disturbed  by  the  half- sine  pulse  of  white 

noise.  Such  responses  were  continuously  monitored  and  the  positive  and 

negative  peak  to  rms  values  were  recorded  and  used  to  form  P  (3)  versus  3 

M 

plots  in  a  manner  identical  to  that  for  the  rectangular  pulse. 

The  results  of  the  analog  study  are  summarized  in  Figures  15  through 
20.  These  figures  are  families  of  curves  in  P  (P)  for  both  stationary  and 
nonstationary  conditions  with  the  peak  to  rms  response  ratio  (P)  plotted  versus 
the  dimensionless  time  parameter  (T  ).  Figures  15,  16,  and  17  are  the 
results  for  Q  =  5,  20, and  50  respectively,  when  the  nonstationary  excitation  is 
the  rectangular  pulse  of  white  noise.  Figures  18,  19,  and  20  are  the  results 
for  Q=  5,  20,  and  50  respectively,  when  the  nonstationary  excitation  is  the 
half- sine  pulse  of  white  noise. 

All  of  the  curves  exhibit  the  same  characteristic  behavior:  that  is  a 
very  rapid  rise  in  |3  for  the  initial  intervals  of  T  ‘  ,  passing  through  a  knee 

s'c 

in  the  vicinity  of  T  =  1,  then  increasing  in  P  very  slowly  beyond  the  knee. 

As  would  be  expected,  the  magnitudes  of  the  curves  for  the  pulsed  nonstationary 

excitation  are  bounded  by  the  stationary  curves  from  Section  3.2.  These 

figures,  therefore,  can  be  used  to  assess  the  time  duration  for  an  oscillator 

* 

to  effectively  achieve  stationarity  in  its  response.  Beyond  T*  =  1,  the  p 
results  for  stationary  white  noise  can  be  used  for  predicting  peak  to  rms 
response  values  and  are  noted  to  yield  conservative  estimates  for  p  in  all 
cases . 
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Figure  13.  Typical  Excitation  of  the  Mechanical  Oscillator  Subjected 
to  Pulsed  Random  Excitation;  Half-Sine  Envelope 


Figure  14.  Typical  Responses  of  the  Mechanical  Oscillator  Subjected 
to  Pulsed  Random  Excitation;  Half-Sine  Envelope 
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Figure  15.  Ratio  of  Peak  Response  to  RMS  Response  of  the  Mechanical  Oscillator  (p)  Versus  the  Dimensionless 
Time  Parameter  (T^)for  Stationary  and  Pulsed  Random  Excitation.  Rectangular  Envelope,  Q  =  5 


Ratio  of  Peak  Response  to  RMS  Response  of  the  Mechanical  Oscillator  (p) versus  the  Dimensionless 
Time  Parameter  (T*)  for  Stationary  and  Pulsed  Random  Excitation.  Rectangular  Envelope,  Q  =  20 
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Figure  17.  Ratio  of  Peak  Response  to  RMS  Response  of  the  Mechanical  Oscillator  (p)  versus  the  Dimensionl 
Time  Parameter  (T#)for  Stationary  and  Pulsed  Random  Excitation.  Rectangular  Envelope,  Q  =  5 


Figure  18.  Ratio  of  Peak  Response  to  RMS  Response  of  the  Mechanical  Oscillator  ((3)  versus  the  Dimensionless 
Time  Parameter  (T*)  for  Stationary  and  Pulsed  Random  Excitation.  Half-Sine  Envelope,  Q  =  5 


Figure  19.  Ratio  of  Peak  Response  to  RMS  Response  of  the  Mechanical  Oscillator  ((3)  versus  the  Dimensionless 
Time  Parameter  (T*)  for  Stationary  and  Pulsed  Random  Excitation.  Half-Sine  Envelope,  Q  =  20 


20.  Ratio  of  Peak  Response  to  RMS  Response  for  the  Mechanical  Oscillator  (|3)  Versus  the  Dimensio 
Time  Parameter  (T*)  for  Stationary  and  Pulsed  Random  Excitation.  Half- Sine  Envelope,  Q  =  50 


4.  CONCLUDING  REMARKS 


4.  1  SUMMARY 

This  report  considers  the  peak  response  behavior  of  a  linear  mechani¬ 
cal  oscillator  when  subjected  to  stationary  and  nonstationary  random  excitation. 
The  peak  response  is  expressed  as  a  ratio  of  peak  to  rms  response  and  is 
denoted  as  (3.  The  stationary  random  excitation  is  white  noise  and  the  non¬ 
stationary  random  excitations  are  pulses  of  white  noise  shaped  in  time  as 
either  a  rectangle  or  a  half-sine.  The  peak  response  results  were  obtained 
by  an  analog  computer  study  and  are  displayed  as  a  family  of  constant  proba¬ 
bility  curves  with  (3  plotted  versus  a  dimensionless  time  parameter  denoted 
❖ 

as  T  .  This  parameter  contains  the  undamped  natural  frequency  of  the 
oscillator  (f^),  the  damping  of  the  oscillator  expressed  as  Q,  and  a  time 
duration  noted  as  T  for  stationary  white  noise  and  At  for  the  shaped  white 
noise. 

The  (3  results  for  stationary  white  noise  are  seen  to  establish  an  upper 
bound  on  the  (3  results  for  the  pulses  of  shaped  white  noise.  Thus,  the 
stationary  white  noise  (3  plots  can  be  applied  to  conservatively  predict  the 
peak  response  of  an  oscillator  to  pulsed  random  excitation.  For  T  values 
<1  (that  is,  the  time  duration  of  the  pulse  is  relatively  short  as  compared 
to  the  natural  period  of  the  oscillator),  the  (3  results  for  stationary  white 
noise  yield  highly  conservative  estimates.  For  T  ^  1,  the  white  noise  (3 

results,  although  still  conservative,  are  more  nearly  the  same  order  of 
magnitude  as  the  (3  results  for  the  shaped  white  noise.  For  a  small  number 
of  response  cycles,  the  peak  to  rms  response  ratio  is  noted  to  be  dependent 
on  the  Q  of  the  oscillator  for  each  For  a  larger  number  of  response 

cycles,  however,  these  peak  to  rms  results  tend  to  become  independent  of  Q. 

F or  the  linear  oscillator,  the  f3  curves  provide  an  empirical  answer 
to  the  following  basic  nontrivial  questions  which  are  of  practical  importance 
in  mechanical  design: 
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(1)  How  high  an  amplitude  would  be  reached  in  a  finite  time  inter¬ 
val  T  with  probability  ^  is  recalled  that  P^.(|3)  is 

defined  as  the  probability  that  M,  the  maximum  value  for  the 
envelope  peak  in  time  T,  is  <  (3  times  the  rms  response  value. 

(Z)  What  length  of  time  is  necessary  for  an  oscillator  to  effectively 
achieve  stationarity  in  its  response  after  being  excited  by  a 
pulse  of  random  excitation? 

The  first  question  arises  in  extreme  value  problems  where  one  is  concerned 
with  the  likelihood  of  a  single  catastrophic  event.  For  example,  what  is  the 
probability  that  some  basic  structure  of  an  aircraft  will  experience  an  extreme 
stress  beyond  ultimate  strength  due  to  a  gust  load  during  its  service  life? 

For  a  second  example,  what  is  the  probability  that  a  resiliently  mounted  equip¬ 
ment  package  in  a  spacecraft  will  collide  with  neighboring  structure  due  to 
random  vibration  during  the  launch  phase?  The  second  question  arises  in 
those  cases  where  one  is  concerned  with  making  stationary  approximations 
for  nonstationary  environments.  For  example,  will  a  stationary  vibration 
test  properly  simulate  the  most  severe  response  conditions  which  a  com¬ 
ponent  experiences  during  exposure  to  a  nonstationary  vibration  environment. 

Continuing  with  the  example  of  vibration  testing,  assume  the  first 
crossing  of  a  given  extreme  stress  level  will  produce  a  failure  of  the  component 
to  be  tested.  Such  a  failure  criterion  is  often  acceptable  for  components  sub¬ 
jected  to  intense,  short  duration  vibration  environments  such  as  launch  environ¬ 
ments  for  spacecraft.  As  indicated  by  the  (3  curves,  the  probability  of  a 
given  extreme  value  is  a  function  of  both  the  rms  value  for  the  vibration  and 
the  exposure  time.  Hence,  the  (3  curves  can  be  used  as  a  criterion  to 
arrive  at  appropriate  test  levels  by  specifying  (1)  the  time  varying  rms 
magnitude  of  the  environment,  (2)  an  acceptable  peak  value  probability, 
i.  e. ,  a  P^.(|3)  value,  and  (3)  the  expected  time  exposure  in  the  random 
environment.  From  these  data  and  the  (3  curves,  a  stationary  vibration 
test  can  be  specified  which  simulates  a  nonstationary  vibration  environment. 
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The  (3  curves  depict  in  a  general  way  the  time -probabilistic  be¬ 
havior  of  the  peak  response  of  a  linear  oscillator  excited  by  both  stationary 
and  nonstationary  random  excitation.  These  results,  however,  can  be 
applied  to  provide  qualified  solutions  to  various  specific  problem  areas  such 
as  for  test  specifications  which  are  discussed  briefly  on  the  previous  page. 

It  must  be  understood  that  conclusions  drawn  from  the  (3  curves  are  based 
upon  an  extreme  value  criterion  and  can  be  applied  to  problems  only  where 
such  a  criterion  is  plausible. 

4.  2  RECOMMENDATIONS  FOR  FUTURE  STUDIES 

On  the  basis  of  this  report,  several  important  topic  areas  are 
recommended  for  future  study.  For  the  mechanical  oscillator: 

•  Repeat  the  analog  study  for  other  values  of  f  ,  Q,  and  T 

n 

(or  At)  to  establish  response  results  which  are  statistically 
more  significant. 

•  Consider  envelope  shape  functions  other  than  a  rectangle  or 
half- sine  (for  example,  a  triangle  or  a  trapezoid). 

•  Consider  the  effects  on  the  (3  parameter  of  varying  the 
oscillator  center  frequency  (f  )  over  the  frequency  band  of 
the  input  excitation  where  the  input  spectrum  is  other  than 
white  noise. 

•  Consider  the  effects  of  nonlinearities  on  the  (3  parameter. 

Similar  studies  also  should  be  performed  for  multi  degree-of-freedom 

systems  to  determine  the  significant  parameters  affecting  the  (3  response. 

Although  the  above  recommendations  can  be  effectively  carried  out 

using  electrical  analog  equipment  (as  in  this  report)  an  alternative  program 

would  be  to  perform  the  suggested  analog  experiments  by  digital  or  hybrid 

methods.  In  these  ways, the  vast  amount  of  peak  response  data  can  be 

routinely  and  efficiently  organized  to  form  the  desired  P  (|3)  versus  (3  plots 

as  well  as  the  (3  versus  f  T/Q  curves. 

n 
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